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ABSTRACT 

We present Spitzer observations of a sample of 12 starless cores selected to have prominent 24 /jm 
shadows. The Spitzer images show 8 and 24 /urn shadows and in some cases 70 /im shadows; these 
spatially resolved absorption features trace the densest regions of the cores. We have carried out a 
12 CO (2-1) and 13 CO (2-1) mapping survey of these cores with the Heinrich Hertz Telescope (HHT). 
We use the shadow features to derive optical depth maps. We derive molecular masses for the cores 
and the surrounding environment; we find that the 24 /im shadow masses are always greater than or 
equal to the molecular masses derived in the same region, a discrepancy likely caused by CO freeze-out 
onto dust grains. We combine this sample with two additional cores that we studied previously to 
bring the total sample to 14 cores. Using a simple Jeans mass criterion we find that ~2/3 of the cores 
selected to have prominent 24 /im shadows are collapsing or near collapse, a result that is supported 
by millimeter line observations. Of this subset at least half have indications of 70 pan. shadows. All 
cores observed to produce absorption features at 70 /im are close to collapse. We conclude that 24 /im 
shadows, and even more so the 70 /im ones, are useful markers of cloud cores that are approaching 
collapse. 

Subject headings: Dust, Extinction, ISM: Molecules, ISM: Globules, ISM: Clouds, Stars: Formation 



1. INTRODUCTION 

Stars are b orn in cold cloud cores (e.g., 

iBergin fc Tafallal I2007T ) . where gas and dust are 
compressed to high enough densities to cause the core 
to start collapsing. One of the most critical questions 
in astronomy is what the physical conditions for core 
formation and collapse are. Therefore, observations 
of cores in the rare stage during, or just prior to, 
collapse are critical to constrain theo ries of the very 
earliest stages of star-for matio n (e.g.. IShu et al.l 119871 : 
iBallesteros-Paredes et aHl200a lMversll2005l ). 

The initial process of collapse, or the transition from a 
stable core to a core with a n embedded protostar, hap- 
pens rapidly (|Havashilll966l ): additionally, the cores are 
necessarily dense and cold. These conditions pose sig- 
nificant observational challenges to identify cores close 
to collapse. Because of the temperatu r e ranges in- 
volved, ~ 10 K (e.g.. iLemme et afl 119961 : ICaselli et al.1 
ll999t lHotzel et al.ll2q02D. a nd high densities. ~10 5 cm" 3 
(e.g-jBacman n et al.ll2000D . cold cloud cores can best be 
observed at far-in frared, submillim eter, and millimeter 
wavelengths (e.g-. IStutz et al.ll2009f ). Such observations 
of dense cores show that most o f them are close to equi- 
librium and not collapsing (e.g.. lLada et al. 2008). Fur- 
thermore, if the cores are not supported by thermal and 
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turbulent pressure alo ne, then a modest magnetic field 
can halt collapse (e.g., iKandori et al.l 120051 : IStutz et all 
[2007tlAlves etalllBol . 

Millimeter line observations are the traditional way to 
search for collapsing cores (e. g. . 1 Walkerlt7aT]fl986T l : how- 
ever, there are ambigu ities in the interpr e tation of such 
meas urements (e.g . . IWalker et al.l 119881: iMenten et al.1 
Il987t iMundv et aLlll990t iNaravanan et al.lll998l ). Here 
we use mid-infrared shadows as an alternative observa- 
tional approach to the study of starless cores, one that 
does not depend on the interpretation of millimeter line 
profiles and how these profiles relate to the underlying ve- 
locity field of the core material. Additionally, our method 
is sensitive to higher column densities than traditional 
studies of starless cores by background star extinction, 
usually limited to Ay < 30 magnitudes (e.g.. lAlves et al.l 
[20011: IKandori et al.ll2005h . 

Cores that are isolated from regions of massive star for- 
mation are very useful to understand the initial stages of 
collapse into stars. These cores are free from complicat- 
ing factors, e.g., massive core fragmentation, feed-back 
from high- mass st a rs, and protostellar outflows (e.g., 
Lada fcLadal 120031 : IdeWit et al.l 120051 : iFallscheer efatl 



2009), to name a few, that can have large effects on a 



study aimed at understanding how individual low mass 
stars form. We present a sample of 12 relatively isolated 
cores that were selected from Spitzer MIPS observations 
of Bok globules and other star-forming regions to have 
prominent 24 /im shadows, spatially resolved absorption 
features caused by the dense core material viewed in ab- 
sorption against the interstellar radiation fie ld. We also 
inclu de 2 previou sly studied cores, CB190 (|Stutz et all 
I2007h and L429 (jStutz et all 120091 ). for a final sample 
of 14 cores. Such objects with 24 /im shadows always 
show a counterpart 8 /im shadow, analogous to the IR 
absorption features produced by more distant and mas- 
sive structures termed infrared dark clouds (IRDCs; e.g., 
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Perault et aT1ll996HButler & Tan|[200l iPeretto fc Fulled 
20091: iLee et alJl2009h and the IS O 7 urn absorption fea- 
tures studied bv lBacmann et all (|2000l ). We also present 
Heinrich Hertz Telescope (HHT) 12 CO (2-1) and 13 CO 
(2-1) on-the-fly (OTF) maps of regions ~ 10' on a side 
surrounding each core. 

Our search method is efficient at identifying cores that 
are near collapse. In § 2 we describe the observations 
and data processing; in § 3 we derive optical depth maps 
from the 8 and 24 /jm images; in § 4 we present our mass 
measurements and we describe the stability analysis that 
we apply to the sample of cores; and finally in §5 we 
summarize our main conclusions. 

2. OBSERVATIONS AND PROCESSING 

2.1. Spitzer IRAC and MIPS data 

The 12 cores presented here were selected from two 
observing programs: Spitzer GTO program ID 53 and 
GO program ID 30384. Core names, MIPS and IRAC 
AOR numbers, and corresponding program IDs are listed 
in Table 1. The IRAC data consist of observations in 
all four channels from programs 139 (P.I. N. Evans), 
94 (P.I. C. Lawrence), and 20386 (P.I. P. Myers). The 
frames were processed using the Spitzer Science Center 
(SSC) IRAC Pipeline vl4.0, and mosaics were created 
from the basic calibrated data (BCD) frames usin g a 
custom IDL program; see iGutermuth et alJ (2008) for 
details. The MIPS observations at 24, 70, and 160 (im 
are from programs 53 (P.I. G. Rieke) and 30384 (P.I. T. 
Bourke). All observations were carried out in scan map 
mode and were reduced u sing the Data Analysis Tool 
(DAT: iGordon et all [20051) . These data w ere pr ocessed 
according to steps outlined in lStutz et al.1 (|2007f ): please 
refer to that publication for further details. The data for 
each core are shown in Figures [T] through Q21 

2.2. 12 CO (2-1) and 13 CO (2-1) Data 

These cores were mapped in the J = 2-1 transitions of 
12 CO (230.538 GHz) and 13 CO (220.399 GHz) with the 
10-m diameter Heinrich Hertz Telescope (HHT) on Mt. 
Graham, Arizona. Observations taken in 2007 (CB27 
and CB44) were obtained using a receiver that was a 
dual polarization SIS mixer system operating in double- 
sideband mode with a 4 - 6 GHz IF band. The obser- 
vations taken after 2007 were obtained using a proto- 
type 1.3mm ALMA sideband separating receiver with a 
4 - 6 GHz IF band. The 12 CO (2-1) line was placed 
in the upper sideband and the 13 CO (2-1) one in the 
lower sideband. The spectrometers were filter banks 
with 128 channels of 250 KHz width and separation. 
At the observed frequencies the spectral resolution was 
~0.3 km s _1 , the band- width was ~40.0 km s _1 , and the 
angular resolution of the telescope was 32" (FWHM). 

The cores were mapped with on-the-fly (OTF) scan- 
ning in RA at 10" s _1 , with row spacing of 10" in declina- 
tion. System temperatures were calibra ted by the stan- 
dard ambient temperature load method (|Kutner fc Ulichl 
1981) after every other row of the map grid. Atmo- 
spheric conditions were generally clear and stable, and 
the system temperatures were nearly constant through- 
out the maps. Data for each CO isotopologues were 
processed with the CLASS reduction package (from the 
Observatoire de Grenoble Astrophysics Group), by re- 
moving a linear baseline and convolving the data to a 



square grid with 10" grid spacing (equal to one-third the 
telescope beamwidth). The intensity scales for the two 
polarizations were determined from observations of stan- 
dard sources made just before the OTF maps. The grid- 
ded spectral data cubes were proce ssed with the Miriad 
software package (jSault et al.lll995h for further analysis. 
For each sideband the two polarizations were averaged 
together yielding the map RMS noise levels listed in Ta- 
ble 2; also listed are the map center coordinates, dates 
of observation, sizes, system temperatures, and veloc- 
ity ranges which include significant emission. The 12 CO 
and 13 CO integrated intensity images, integrated over 
the velocity ranges listed in Table 2, are shown in Fig- 
ures Q2] through [15l these images are centered on the 
24 /im shadow coordinates. 

3. ANALYSIS 

3.1. Optical Depth Maps 

In this section we calculate the 24 /im and 8 /im optical 
depths for the shadow features seen in Figures [T1H21 We 
assume that the cores are optically thin and that there is 
no foreground emission originating from the surface lay- 
ers. We app l y an analysis analog ous to that p r esented in 
IStutz et all P007f ) at 24 /im and lStutz etail (|2008bh at 
8 /im; however, we update the analysis presented in that 
work with a more robust calculation of the unobscured 
background flux level (/o in what follows). 

3.1.1. 24 p,m Optical Depth Calculation 

Under the above assumptions, the optical depth at 
24 /im is 



Here /dc is the contribution to the image flux level of all 
foreground sources, such as zodiacal emission and instru- 
mental background; /' and /g are the shadow flux and 
the intrinsic, unobscured flux, respectively. In what fol- 
lows / and /o indicate the DC-subtracted values of the 
shadow fluxes: / = /'-/dc and / = /o-/dc- The con- 
tribution of /dc to the overall image flux level is assumed 
to be constant over image areas large compared to the 
size of the shadow while /o is determined in the region 
immediately adjacent to the shadow. Further details on 
the determination of these two quantities is provided in 
the following discussion. However, before proceeding, we 
note that the shadow depth, or contrast, in the absence 
of foreground emission from the surface of the shadow 
itself, is the result of the depth of the absorption feature 
relative to the quantity (/q-/dc); a stronger background 
allows for robust sensitivity to larger r. Furthermore, an 
over-estimate of /dc Wlu artificially boost the shadow 
contrast and cause an over-estimate of r and all depen- 
dent qu antities. Rather th an following the procedure ap- 
plied in IStutz et al.l (|2009h . where the darkest region of 
the L429 shadow was used to measure foe , the approach 
implemented here takes advantage of dark regions in the 
Spitzer images which are not associated with any ab- 
sorption features. These regions provide an independent 
estimate of the foreground emission; furthermore, when 
shadow pixels are darker than, or have a lower value than, 
the foe level, we exclude them from our analysis because 
they are most likely not sensitive to the total column of 
material in that pixel, i.e., they have become saturated. 



3 



Image DC level (fuc ) determination. To measure /dCi 
we select regions in the 24 /im image with no diffuse 
emission. These regions also have little to no detected 
emission at 70 and 160 /im. As a further check on the 
cleanliness of the background estimate, when possible we 
select the background region to overlap partially with, or 
be adjacent to, 12 CO and 13 CO map areas that are free 
from molecular line emission. Additionally, we avoid the 
areas near the edges of the images due to poor signal-to- 
noise and sparse coverage. The background region box 
size is 5' on a side, corresponding to ~ 59, 000 mosaic pix- 
els, which are subsampled by 2 x 2 relative to the physical 
pixels. We set /be equal to the 0.5 percentile pixel value; 
as discussed above, this low value is conservative because 
a higher value may over-estimate the shadow contrast 
relative to f . The error in /dc> cdc, is determined by 
simulating ~ 59, 000 pixels drawn from a normal distri- 
bution with the same mean and standard deviation as the 
background region; the 0.5 percentile pixel value is stored 
and the process is repeated 10 4 times. The standard de- 
viation of the resulting dis tribution values is our estimate 
of ctdc (|Stutz et al.ll2007l ). The pixel error map for the 
backgroun d subtracted images is then calculated using 
the DAT ()Gordon et alj l2005f) pixel error map, multi- 
plied by a factor of 2 to take into account pixel-to-pixel 
correlations due to the sub-sampling, with ctdc added in 
quadrature. 

\J nob soured flux (fo) determination. We apply itera- 
tive Gaussian fitting to the image pixel value distribu- 
tions to estimate fo, the local, unobscured fl ux near the 
shadow , similar to the method developed bv lStutz et al.l 
(2008a) to reconstruct and analyze the noise properties 
in deep Hubble Space Telescope ACS images. In sum- 
mary, we adjust the upper cutoff of the pixel value dis- 
tribution so that the Gaussian fit is least affected by the 
non-Gaussian tail of the high pixel value distribution. 
This is done by finding the upper pixel value cutoff that 
has the narrowest Gaussian width. This fitting method 
is constructed to minimize the effects of sources or shad- 
ows on the estimate of fo. We proceed by constructing 
a pixel value histogram for all pixels within ro of the 
shadow using a bin size of 0.4 //Jy/arcsec 2 . We then 
clip the distribution by excluding all bins with values 
below the mean pixel value in a 30" radius aperture cen- 
tered on the shadow. The shadow pixel-value distribu- 
tion in this region is indicated as a grey histogram in 
Figures [TBHT^l the mean value of the shadow flux is in- 
dicated with a solid line. Excluding these bins avoids 
a possible bias of the best-fit mean pixel value towards 
lower values. Next we fit the histogram data to an un- 
weighted Gaussian model, using standard IDL packages, 
and excluding all bins with values greater than the maxi- 
mum bin plus one. We repeat this fitting procedure as we 
include successively higher value bins one at a time. The 
resulting best-fit Gaussian width (<j) values are shown 
in Figures [T6H191 in the lower-left panel for each object, 
plotted as a function of the maximum bin value included 
in the fit. The solid line indicates the minimum a from 
the distribution of Gaussian widths; the mean value of 
the fit corresponding to the minimum value of sigma is 
the adopted value of fo ■ These two parameters are listed 
in the figure. 

To determine the effect the choice of pixel histogram 
radius (ro) has on our derived value of fo we repeat this 



procedure over a range of ro values from ro — 60" to 
300", in 20" steps; the resulting values of fo and a are 
shown in the right-hand-side panel of Figures [TfirtTTjl In 
most cases the fo profiles flatten out by ro = 240", as 
indicated by the dashed line. This value is used for all ob- 
jects except CB42, which lies near the edge of the MIPS 
mosaic; since the fo profile for CB42 appears to flatten 
out by 140", the largest value of ro that we test for this 
core, we assume this value is valid. In the cases of CB26 
and CB27 we do not measure as abrupt a flattening. For 
these cores we still use ro = 240" as a conservative esti- 
mate; larger values for ro may over-estimate fo and the 
derived optical depths. Furthermore, ideally we calcu- 
late fo as close as possible to the shadow features; as 
the distance from the core is increased, other effects can 
bias the determination of /o, such as diffuse emission not 
associated with the core. We estimate the uncertainty 
in fo to be er(/o) = 0.02 x f . The Gaussian error in 
the mean gives a SNR for fo of > 230 for all cores; we 
consider this SNR to be optimistic. 

Using our calculated value of fo and the background 
subtracted image, we calculate T24 in each pixel and the 
associated error cr T ,24, using standard propagation of er- 
rors, including ctdc and cr(/o)- The pixel-to-pixel cor- 
relations are significant and we do not take this into ac- 
count in the r pixel error maps; consequently, the error 
maps underestimate the true errors and we do not use 
them in the subsequent analysis. In the cases where the 
values of / near the center of the cores are negative, 
causing T24 to diverge, we set the value of T24 equal to 
the mean (T24) in an aperture with radius 10" or larger, 
depending on the size of divergent region, centered on 
the core. For most (8/12) cores, the number of diver- 
gent pixels is less than 90, corresponding to an effective 
circular area with radius smaller than 6"5; for objects 
CB130 and CB180 the divergent region has an effective 
radius of ~ 7"5, for CB44 the region has a radius of about 
25", and for L1552 the region radius is about 20". This 
approach is conservative when analyzing core stability 
because it effectively underestimates the 24 /im shadow 
masses. Furthermore, an alternative approach is to esti- 
mate fuc as a fraction of the shadow pixel values; while 
this method would circumvent the divergent pixel issue, 
it would yeild larger core masses, potentially biasing our 
collapse analysis. The divergent pixels found with our 
method are flagged and accounted for in the subsequent 
analysis. The 24 /im r pixel maps are shown in the left 
column of Figures Q2] through [T5l 

3.1.2. 8 11m Optical Depth Calculation 

We apply the same techniques described above for the 
24 /itm images to the 8 /xm images, with the exceptions 
described here 7 . The results are illustrated in Figuresl2"0T- 
|2"51 For the /dc determination, it is not possible to 
use the same regions as those used for the 24 /im im- 
ages because of the limited coverage of the 8 /zm images; 
therefore we select new regions that are near the 24 /im 
regions when possible. In all cases the 8 /im /dc re- 
gions are selected to include the darkest regions of the 
image; we set /dc equal to the 0.1 percentile. The back- 
ground region box size is 5' on a side, corresponding to 
~ 60, 000 pixels at the plate scale of the IRAC images 

7 There is no IRAC coverage of source CB42 
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(1'.'22 pix -1 ). The 8 /im error is assumed to be 2% of 
the pixel flux value; the background error is estimated as 
described above. These two quantities added in quadra- 
ture make the background-subtracted pixel error map. 
The method used to determine fo is identical to that de- 
scribed above, the only exception being that for some 
objects (CB23, L1544, L492, and L694-2) the maximum 
radius for which we can measure the pixel value distribu- 
tion is limited by the size of the images. For these four 
objects, we chose the largest radius measurement of fo for 
which over 50% of the area of the annulus contains image 
pixels. For these cores the maximum radius is ro = 240" 
for CB230, and r = 200" for the rest. For cores with 
sufficiently large images, we observe a flattening of the 
fo profiles at a radius of ro = 260"; the adopted values 
for ro are indicated as dashed lines in the right column 
of Figures [HH231 In the cases of cores CB23 and CB27, 
we do not observe a flattening; therefore the t$ maps for 
these sources may be less reliable. 

As described above for the 24 /im images, we use our 
measurements of /q, /', and fnc measured form the 
8 /imimages to calculate 

, ( f — /DC 
t s = -In — 

\Jq~ /DC , 

in each pixel; we also calculate ct Ti s with standard prop- 
agation of errors and including ctdc an( i c(/o), where 
cr(/o) = 0.02 x f . Divergent pixels are assigned a value 
according to the procedure outlined above, are nagged, 
and are accounted for in subsequent analysis. 

We detect 3.6 /Ltm emission above the sky level in 
a large fraction of our cores (TMC1-1C, CB23, CB26, 
CB27, L1544, L1552, and L694-2, and marginally in 
L492, CB130, and CB180). This emission is spatially 
coincident with the shadows, and similar in extent to 
the 160 /im emission. The lack of any significant 8 /jm 
aromatic feature indicates that this emission is likely 
to be caused by the surface layer of dust grains scat- 
tering the interstellar radia tion field, analogous to th e 
near-IR emission studied by iLehtinen fc Mattilal (fl996h . 
Such near-IR e mission has been used to study interstel- 
lar c l ouds (e.g. iFoster fe Goodman] 120061: IPadoan et al.l 
120061: iNakaiimaet al.l I2008D . Uuvela et all (|2008h shoV 
that when using J, H, and Ks imaging this technique 
is sensitive to densities of 1 < Ay < 15 magnitudes. 
Including the 3.6 /im emission might increase the den- 
sities accessible to this technique, as well as providing 
a high-spatial resolution shadow-independent extinction 
estimator. We note that for the cores L1544, CB44, and 
L694-2 the emission shows a limb-brightening geometry 
at 3.6 /im, where the center of the shadow region appears 
darker than the surrounding region. 

The 8 /im shadow determinations provide a test of 
the validity of those at 24 /im. In all cases where a 
good comparison is possible, the agreement is satisfac- 
tory if we take the optical depths at the two wave- 
lengths to be approximately equal, T24/T8 = 1.0 ± 
0.1 (see also iFlaherty et al.1 120071 : rChapman et "all 120091 : 
IChapman fc Mundvl l2009f ) . Additional details will be 
provided in a future publication. 

3.2. Distance measures 

Here we make some remarks on our adopted distances, 
which are summarized in Table 3. 



TMC1-1C, CB23 l CB26, CB27, L1544, L1552 — 

iKenyon et al.l (|1994h derive a distance to the northern 
portion Taurus of 140 ± 10 pc based on extinction of 
background stars. This distance is roughly consistent 
with recent dete rminations using tr igonometric parallax 
with the VLBA (jTorres et al.ll2009h . 
CB42 — We find no distance estimate for this source 
in the literature. In the directi on of CB42, the exten t 
of the local bubble is ~ 100 pc (|Lallement et al.l [20031: 
we adopt a distance of 100 pc, a conservative lower limit 
distance. 

CB44 — Iffllton fc Lahullal (fl995h list the 
iBok fc McCarthy! (Il974h assumed distan ce of 400 pc. 
They also list the iTomita et all (Il979l) es timate of 
300 pc, and the lArquilla fc Goldsmith! (|1986f) distance 
estimate of 400 to 600 pc. We adopt a distance of 
400 pc. 

L492, CB130 — IStraizvs et "all (|2003l ) investigate the 
distance dependence of the interstellar extinction in the 
direction of the Aquila Rift. These authors find that the 
distance is 225 ± 55 pc to the front edge and that the 
cloud has a depth of about 80 pc. Because these cores 
are viewed in absorption they are unlikely to be located 
behind much cloud material; therefore we adopt 225 pc 
as the distance t o these cores. 

CB180 — Iffllton fc Lahullal (|1995h list the 
Bok fc McCarthy (Il974h assumed distance of 400 pc. 
Bok fc McCarthy! (|1974l ) do not elaborate on their 



adopted distances. This object is near the Scutum 
association. We adopt 400 pc for the distance to this 
core. 

L694-2 — iKawamura et al.l (|2001[ ) derive a distance of 
230 ± 30 pc using star counts. 



4. DISCUSSION 

4.1. Interpretation of Collapsing Cores 

The stability of a cloud core is usually deduced 
through observa tions of millimeter-wav e molecular emis- 
sion li nes (e.g., iGreeersen et al.l 119971: iMardones et al.l 
19971 iGregersen fc Evans! l2000t fCrapsi et all I200E' 



Williams et all 120061 : iSohn et all 12007J), but the dy- 



namical state of these lines can be affected by ro- 
tation, turbulence, optical depth, and other effects 
that can hide or mimic the indicators for collapse. 
We use infrared shadows, together with millimeter- 
wave spectroscopy, as an alternative, complementary, 
and reliable ind i cator of impending collapse (e.g., 
IStutz et all 120071 l2009fh These observational tech- 
niques can be used to estimate the masses and dy- 
namical states of the cores. Conventionally, collapse 
is identified by fitting such observations to Bonnor- 
Eber t model profiles which assume spherical symmetry 
(e.g.JEvans et alJ 120011 : iBallesteros-Paredes et all 120031 : 
IShirlev et al.l I2005TK However, as shown in Figures [T] 
through 1 121 few of our objects appear symmetric on the 
plane of the sky, calling the applicability of this approach 
into question. 

Instead, we have probed the collapse state of these 
cores on the basis of a simple Jeans mass calculation. 
In what follows, we determine 24 /im core masses and 
sizes. Given these core properties, and assuming a tem- 
perature of 10 K (reasonable for starless cores) we derive 
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a Jeans mass for each core: 

Mj=(-^-\ m ''r^M- 1 ' 2 , (2) 
V &(J, p m H J 

where M is the 24 /zm core mass (M%£ c in Table 4), 
R is the effective radius of the core (R^ c in Table 4), 
measured from the 24 //m image, and we have assumed 
that the mean molecular mass per free particle \i p is equal 
to 2.4. The sizes of the cores are derived from the area 
in which we detected the 24 /im shadows (see § 4.2.1). 
We note that assuming a temperature of 8 K or 12 K 
for these cores cause a ~30% change in the calculated 
Jeans mass. The distance dependence of the Jeans mass 
is D 1 / 2 . In the following analysis (see Tables 4 and 5) 
we compare shadow masses, with a distance dependence 
of D 2 , to the Jeans mass. We caution that the ratio of 
the 24 /im masses to the Jeans mass has a steep distance 
dependence of D 3 / 2 . This Jeans mass analysis avoids 
possible biases introduced by imposing more stringent 
geometrical assumptions in the core analysis. However, 
it only takes thermal support into account. 

Under more realistic conditions, turbulen ce might also 
provid e support against collapse. From iHotzel et al.l 
( 2002]), the energy contributed by turbulence to the sup- 
port of the cloud is 

Eturb = Jmu - Eth > (3) 

where AMU is the atomic mass number of the species 
from which we measure the observed line width AV, and 
AVth is the expected thermal line width for that species, 
assuming a temperature of 10 K. To estimate possible 
contributions of turbulence to the support in the cores, 
we use previously measured AV values for hi gh density 
trace r s from the literature when available ( e.g. I Lee et all 
l200i iCrapsi et al.ll200fit iSohn et alJ I2007T ) . When such 
measurements are not available, we use other line widths 
as estimators, noting that lower density tracers will over- 
estimate the line width and therefore the contribution 
from turbulence to the support. Often, the turbulent 
support is much less than the thermal support and hence 
can be ignored in considering the stability against col- 
lapse. In the cases where turbulent support is poten- 
tially significant, we attempt to assess its contribution 
relative to the Jeans mass and magnetic field required 
for stability. 

Magnetic fields may also contribute to the support. 
The inferred magnetic fie ld strengths in cores ran ge from 
< 15 yG to ~ 1 60 nG ijBergfn fc Tafallal 12001 . From 
IStahler fe Pallal (|2005l ). given a magnetic field strength 
B and a cloud radius R, the mass that can be supported 
is 

Using this relation, and a measurement of the core 
mass and size, we estimate the minimum magnetic field 
strength required to support the core. If the cores re- 
quire field strengths above ~ 160 /iG, we consider it 
unlikely that the core is supported by a magnetic field. 
Furthermore, if the minimum magnetic field derived in 
this fashion is reasonable, but the core mass is signifi- 
cantly greater than the Jeans mass (by a factor > 2), we 



consider it unlikely that the core is stable against col- 
lapse. Our understanding of the magnetic field strengths 
in cores does not allow for strong constraints. However, 
this qualitative assessment gives a useful estimate of the 
possible role of the magnetic field. 

Finally, rotation can also contribute to the support in 
cores. In column 9 of Table 4 we list the rotational 
velocity required to fully support our sample of cores. 
Future high-resolution observations will determine the 
extent of the role of rotation in these cores. In the anal- 
ysis that follows, we interpret the line widths as caused 
only by turbulent and thermal motions; more realisti- 
cally, turbulen t motions on large scales can b e confused 
with rotation (|Burkert fc Bodenheimerl l2b00). Further- 
more, irregular infall, or fragments, can mimic turbu- 
lence; inspection of Figures fTI through [T2l reveals that a 
significant fraction of cores have geometries which are not 
regu lar, reminiscent of structure s seen in s elf-gravitating 
2D (iBurkert fc Hartmannl I2004D and 3D (|Heitsch et all 
2008) simulations. Because of these irregular shapes, we 
have attempted to minimize geometrical assumptions in 
the analysis of our sample of cores, as discussed above; 
however, the stability arguments applied here inherently 
still assume approximate sphericity. 

4.2. Mass Measurements 

We describe the method used to calculate masses in the 
cores using the 24 /jm shadows, and in the regions sur- 
rounding the cores using the 12 CO and 13 CO data. The 
results are summarized in Table 4. We also list the Jean 
mass for each core, assuming all the cores are at a tem- 
perature of 10 K and that the core density is described by 
the shadow mass (M^f 6 ) and the effective radius (R^ ; 
see below). 

4.2.1. Core Masses: Shadow Masses 

We calculate mas ses in the c o res ap plying a similar 
method as that of IStutz et "all (|2007| ). Core shadow 
masses, assumed distances, and effective areas are listed 
in Table 4. In these calculations, we use the 24 11m im- 
ages a nd 8 /im dust opacities from lOssenkopf fc Henninj 
(1994): k 8 = 1.21 x 10 3 cm 2 gm" 1 for thick ice mantles 
at a density of 10 6 cm~ 3 . This choice of dust opac- 
ity is motivated by our analysis in § 3.1.2 where we 
find that r 2 4/g = T2a /t% = 1.0 ± 0.1. This choice of 
Ossenkopf fc Henningl (|1994h dust model has a ratio of 
T24/T8 = 1.1, close to our measured value. 

The mass in a region with an optical depth T24 is: 

flD 2 

M 24 = a x r 24 , (5) 

where is the solid angle subtended by the region, D 
is the distance, a(= 100) is the assumed gas-to-dust 
ratio, and k is the assumed dust opacity. We measure 
the mean flux in an aperture of 1.'5 radius, excluding all 
pixels that have T24 < 0.0 in the maps described above. 
Using the measurements of /be and /o and the mean 
flux per pixel in the shadow, we then calculate the mean 
T24 per pixel. Using equation 2 and the number of pixels 
(fi = N p i X x pixel scale 2 ) included in the mean flux cal- 
culation we calculate a shadow mass. We use the value 
of Q in subsequent analysis as a measure of the effec- 
tive core sizes, where R c °ff = The uncertainties 
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in the masses measured in this fashion are of order 15%, 
without including distance uncertainties. We cross-check 
these masses by calculating core masses directly from the 
r-pixel maps. We find that the two sets of masses agree 
well; however, the mass errors are extremely low for the 
r-pixel map masses. We attribute this discrepancy to 
highly correlated errors, as discussed in § 3.1, and disre- 
gard the small uncertainties in favor of the ~ 15% errors. 
The agreement between the two masses indicates that 
the method we use to estimate foe (see § 3.1.1) does 
not have a large effect on the calculated core masses, 
even accounting for divergent pixel values. 

4.2.2. Envelope Masses: Molecular Masses 

We use the 12 CO (2-1) and 13 CO (2-1) maps to de- 
rive the total column density of hydrogen, N(H) in the 
area surrounding the cores. As listed in Table 2 and 
shown in Figures [131 through II 5 | the CO maps cover ar- 
eas that are much larger than the 24 /xm shadow area. 
The effect of CO freeze-out onto d ust g rains is severe 
(jBacmann et all 120021 : iTafalla et all 12002); therefore, in 
the area where we observe 24 /im shadows, the column 
densities derived from CO observations are likely to be 
underestimated. However, in the regions outside of the 
shadows, where the densities are much lower than in the 
cores themselves, the CO-derived N(H) column should 
have a much higher degree of fidelity to the true column 
in the observed region. Therefore, we use the CO maps 
to derive masses of the larger environment surrounding 
the cores, observed as 24 /im shadows 8 . 

To derive m asses from the CO d ata, we foll ow a similar 
proce dure as iPovich et al.1 (|2009f ) (see also iKang et all 
2009) and we apply the es cape probabi li ty ra diative 
transfer model for CO from iKulesa et al.l (|2005) . The 
12 CO (2-1) and 13 CO (2-1) data cubes are used to de- 
rive the column of CO in each 10" x 10" map pixel, 
assuming an isotope ratio of 12 CO/ 13 CO = 50. The 
model uses the 12 CO and 13 CO observations to solve for 
the column density of CO, N(CO), over a model grid of 
temperatures, densities, and column densities; given an 
initial guess for the density this procedure then fits si- 
multaneously for the temperature and column density in 
eac h map pixel. The mass calcu lation also incorporates 
the lvan Dishoeck fc Blackl (|1988f ) PDR photodissociation 
model of CO, intended to account for the effect of the 
far-UV field, Go, incident on the cores. The final CO 
column that we calculate therefore should include the ef- 
fect of UV dissociation; if the UV-field is not strong, the 
application of this model can cause an over-estimate of 
the column densities. 

After calculating the total (atomic plus molecular) col- 
umn N(H) in each pixel, we derive masses using the dis- 
tances collected from the literature (listed in Table 3 and 
discussed in § 3.2). We calculate total cloud masses by 
integrating the N{H) pixel distribution over a 4' radius 
aperture centered on the shadow, indicated as Mco,4 m 
Table 4. The left column of Figures l24l through l27l show 
the 13 CO integrated intensity maps (grey scale) with the 

8 The 13 CO emission sometimes appears to be bright along an 
edge of the shadow region (seen most clearly in L1552 and CB44) 
and in one case we clearly observe a deficit of emission in the re- 
gion coincident with the shadow (L1544). These characteristics are 
likely caused by the freeze— out of CO on to dust grains, as well as 
by a non-uniform density and temperature structure in the cores. 



calculated N(H) column density contours overlaid; we 
find good spatial agreement between the two. The right 
column of Figures [M] through [27] shows the 12 CO and 
13 CO spectrum averaged over the r < l'.b shadow region 
and in the 4' > r > 1.'5 annulus where the dotted lines 
indicate the velocity range used to calculate the column 
density N(H). We also calculate molecular core masses 
over a circular aperture of size equal to the effective area 
of the core derived from the 24 /im r maps (M™q c ); the 
effective radii (R^ ro ) range from 1.'2 to 1 f 4 for this sam- 
ple of cores and are listed in Table 4. 

4.3. Results 

We use our optical depth maps at 8 and 24 /im, t& 
and T24, to derive core masses; we use the 12 CO and 
13 CO maps to measure the mass of material in the larger 
environment surrounding the cores. We then compare 
the Jeans mass criterion and other information to assess 
the collapse status of each core. Here we summarize our 
results: 

TMC1-1C: Our analysis shows that this core is stable, 
assuming only thermal support. It is therefore highly 
unlikely that this core will collapse soon given its cur- 
rent configuration. If the mean density in the core in- 
creases by ~ 20% the core will exceed the Jeans mass 
criterion; if that happens and if magnetic and turbu- 
lent support are negligible then even tually collapse will 
take place. iKauffmann et al.l (|2008l ) use the MAMBO 
bolometer array to observe the thermal dust contin- 
uum emission at 1.2 mm. In their study this shadow 
is designated as TMC1-1C 3C, with a measured mass 
of 10.71 ± 0.07 M Q , in good agreement w i th ou r CO- 
measured mass of 11.4 M Q . iCrapsi et al.l ((20051) mea- 
sure the N 2 H+ (2-1) and N 2 D+ (2-1) line widths to 
be 0.21 and 0.13 km s _1 , respectively; they do not 
measure a large N 2 D + /N 2 H + ratio and do not classify 
this source a s very advanced chemically or in evolution. 
ICaselli et afl (|2008l ) measure the ortho-H 2 D+ (li )0 -lx,i) 
emission, which is likely correlated with an advanced 
chemical age, and also find that TMC1-1C does not ap- 
pear to be chemically evolved. These results are consis- 
tent with our conclusion that the core is probably stable. 
CB23: Our analysis shows that the 24 /im shadow mass 
exceeds the Jeans mass by ~ 25%, indicating that this 
core is only marginally unstable, if only thermal support 
is considered. The measured N 2 H + line widths (see be- 
low) are narrow and indicate that the turbul ent contribu- 
tion to the support is only of order ~ 30% (|Hotzel et al.l 
120021 ). A magnetic field of - 120 fiG would be suffi- 
cient to support the 24 /xm mass and to halt collapse 
(|Stahler fc Pallall2005| ). A magnetic field of this level is 
on the high side of the range of field stren gths measured 
for other cores dBergin fc Tafallal 120071 ). but remains 
plausible. iPark et al l (| 2004l) observe 12 CO and 13 CO 
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emissi on; the line profiles appear Gaussian. 
(|200l observe CS (3-2) and DCO+ (2-1); 
(1999) observe CS (2-1) and N 2 H+ (1-0). A blue asym- 
metry is observed in CS (2-1); a small amount of a 
blue asymmetry is seen i n CS (3-2 ) and DCO + (2-1) 
although it falls below the iLee et al.l (120041) threshold to 
be classified as a collapse candidate. ICrapsi et al.l (|2005l ) 
measure the N 2 H+ (1-0) and N 2 D+ (2-1) line widths to 
be 0.17 and 0.23 km s , respectively but do not clas- 
sify this source as chemically or kinematically evolved. 
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ISohn et "all (|2007f ) observe an infall asymmetry in the 
HCN (1-0) line profiles of at least 2 of the 3 hyper-fine 
components, and p ossibly in all three. We note that the 
ICrapsi et all (|2005l ) and the ISohn et alj (|2007h observa- 
tions were taken at locations separated by ^60", how- 
ever, both locations are well within the shadow feature. 
Taken with our results, these observations indicate that 
this core is approaching collapse. 

CB26: Our analysis indicates that this core has a mass 
that exceeds the Jeans criterion by a factor of ~3. The 
field required to support the mass of CB26 is of order 
20 /iG, much larger th an the upper limit of 74 /iG set 
by lHenning et all (|2001l ) in a nearby region of CB26. The 
HCN (1-0) line wid ths of ~ 0.6-0.9 km s" 1 measured by 
iTurner et al.l (|1997[ ) indicate that the turbulent support 
is of order ~ 2 to 4 times the thermal support; however, 
ITurner et all (|1997t ) do not show their HCN measure- 
ment and our CO spectra show 2 velocity components, 
which co uld broaden the line widths (see Figure |24|) . We 
note that ITurner et all (|1997l ) do not list the coordinates 
of these observa tions and we assume they ar e the same 
as those listed in lClemens fc Barvainis! (|1988h . which are 
about 1' from the center of the shadow. Without higher 
density tracer line profiles from which to infer the ve- 
locity field of the material, we consider it likely that 
CB26 is near collapse. The young T-Tauri star asso- 
ciated with this source has received far more attention 
than the nearby s tarless core. This YSO has been stud- 
ied in detail (e.g., [L aunhardt et al. 1998; Hcnning 'et alJ 
20011: IStecklum et all 120041 : lLaunhardt et alJ 120091 : 
Sauter et al.l 12009). It lies approximately 3'.5 to the 
south-west of the shadow center, and can be seen in 
Figure [3] as the brightest 24 fim source. Although it 
is not clear what evolutionary, chemical, or other, re- 
lationship the shadow and the YSO might share, it is 
interesting that iHenning et~ al. (2001) measure relatively 
well-ordered submillimeter dust polarization in a region 
~ 40" on a side towards the young star, perhaps indi- 
cating that the magnetic field in the larger region, and 
therefore in the shadow, might also be well behaved. The 
polarization vectors point roughly north-south, indicat- 
ing that the magnetic field is oriented in the east-west 
direction, the same direction as the longer axis of the 
entire CB26 region. 

CB27: Our analysis indicates that this source is stable, 
if supported only by thermal pressure. With an addi- 
tional plausible level of turbulent and magnetic support, 
the core remains stable. This conclusion is supported by 
the molecular l ine observations, which show no evidence 
of asymmetry. iLee et a l. (2004) measure CS (3-2) and 
DCO + (2- 1) line profil e s that are indistinguishable from 
Gaussian. iPark et all (l200l observe 12 C O and 13 CO 
Gaussian profiles as well. ISohn et alJ (|2007l ) observe this 
source to ha ve a red asymmetry in all three HCN (1-0) 
line profiles. iGupta et all (|2009h dete ct C 6 H~ emission 
with a line width of ~ 0.15 km s _1 ; iRoberts fc Millar] 
(2007|) detect N 2 H+ and N 2 D+ emission whose line pro- 
files appear symmetric. 

L1544: We find that our 24 /xm mass is lower than 
the Jeans criterion for collapse; however, the 70 /zm im- 
age indicates that this mass estimate is a lower limit. 
As can be seen in Figure El in the 70 /im image 
there is a hint of a dark region at the location of the 
24 /im shadow. The quality of the 70 /zm image does 



not a llow for a robust analysis; however, IStutz et al.l 
(2009) show that in the case of starless core L429, the 
24 /im shadow mass is underestimated compared to the 
70 yum shadow by a factor ~ 5. Taking this result, 
a plausible corrected 24 /im mass is 7 M Q , a factor 
of ~ 2.5 larger than the Jeans mass, making the core 
much more likely to be approaching col lapse. Many au- 
thors have observed this source (e.g., ILee et al.l 119991: 
Bacmann et alJfeOOCtlShirlev et al.ll2000t]Bacmann et al. l 
20021 iTafalla et alj|2002t IBacmann et alj|2003t ILee et al l 
20031 12004 IWilliams et all l2006t IPark et ; all 12004 
Caselli et al , 120081: iGirota et alJl200flk the consensus is 
that thi s starless c o re is y ery evolved and likely to be col- 
lapsing. ILee et all (|2004l ) observe a CS (3-2) profile with 
a double peak with a brighter blue peak, and a DCO + 
(2-1) profile with a blue asymmetry. ICrapsi et al.l (|2005l ) 
measure the N 2 H+ (2-1) and N 2 D+ (2-1) line widths to 
be 0.32 and 0.29 km s~\ respectively; the N 2 H+ (3-2) 
and N 2 D + (3-2) were 0.19 and 0.34 km s" 1 , respectively. 
They also find that this core has a large N 2 D + /N 2 H + 
ratio equal to 0.23 and classify it as a highly evolved 
core, based on various observatio nal chemical evolution- 
ary probes. IWilliams et al.l (|2006l ) compare the N 2 H+(1- 
0) emission from L1544 and L694-2 and conclude that 
both cores will form stars in ~ 10 4 yr. ISohn et"ai] (|2007f ) 
measure all three HCN (1-0) components to have two 
very strong peaks, all with the blu e p eak higher than 
the re d peak. IBacmann et alJ (|2002f ) and [Bacman n et ahl 
(l2003h find high de grees of CO depletion and a large deu- 
terium abundance; Caselli e t al.l (|2008f ) measure a large 
amount of ortho-H 2 D + (li,o _ emission, one of the 
highest in their sample, and conclude that emission from 
this species is correlated with CO depletion, deuterium 
fractionation, and an advanced chemical evolutionary 
stage. These results are consistent with our conclusion 
that this core is likely to be approaching collapse. 
L1552: In this case, our analysis shows that the 24 /im 
shadow mass exceeds the Jeans mass by a wide margin, 
a factor of ~8. This source also shows a clear detection 
of a 70 /im shadow; applying a correction of a factor of 5 
to th e 24 /zm mass (see LI 544 discussio n andlStutz et ahl 
2009) makes this core very unstable. ILee et alX ( 2004 s ) 
measure a blue asymmetry in both CS (3-2) and DCO + 
(2-1). The ISohn et all lf2007h HCN (1-0) profiles show a 
mild blue asymmetry in the line profiles and report an 
N 2 H + narrow line width of 0.22 km s _1 , indicating that 
turbulent support is negligible. These line measurements 
are consistent with our conclusion that this core is near 
collapse. 

CB42: Our analysis indicates that this source is very 
stable against collapse, with a Jea ns mass that exceeds 
the shadow mass by a factor ~ 10. ITurner et all (|1997f ) 
measure an HCN line width of ^0.82 km s _1 . 
CB44: Our analysis indicates that this core has a mass 
that exceeds the Jeans criterion by a factor of ~ 12. 
Molecul ar line observations of this source are sparse; 
however. lLaunhardt et all (|1998l ) measure a CS (2-1) line 
width of 1.5 km s _1 , with a profile showing some indi- 
cation that it may have a red asymmetry. Observations 
of other (high density) tracers are needed. The magnetic 
field required to support the 24 /im shadow mass is of 
order 190 /iG, large compare d to observations of other 
cores (B ergin fc Ta falla 2007). Based on this evidence 
and without line observations of higher density tracers 



we consider it likely that CB44 is approaching collapse. 
L492: Our mass measurement indicates that this core 
may be unstable to collapse. Although the 24 /xm mass 
only marginally exceeds the Jeans mass, by ~ 15%, 
this source has a 70 /xm shadow, which, as discussed 
above, indicates t hat the 24 /xm mass is a l ower limit 
(jStutz et al.ll2009f) . IKauffmann et all (|2008l ) report a 
1.2 mm dust mass of 18.35 ± 0.17 M Q , in reasonable 
agreement with our CO-measured mass of 23.6 M Q , 
accounting for our larger CO aperture (4'). At both 
8 /im and 24 /im we observe two sub-components of this 
core, oriented in the east- west direction; thi s structure 
is not evident in th e IKauffmann et alJ (|2008f ) MAMBO 
maps. ICrapsi et all {2005) measure the N 2 H + (2-l) and 
N 2 D+(2-l) line widths to be 0.26 and 0.22 km s -1 , re- 
spectively, at a position centered on the western sub- 
clump. The narrow line widths indicate that turbulence 
does not play a significant role in the support of this 
core. They classify this source as marginally evolved be- 
cause it exhibits onl y 3 out of t he 8 e volutionary probes 
that they consider. iLee et all (|2004D measure a CS (3- 
2) with a blue asymmetry, DCO + (2-1) profile with a 
blue asymmetry and possible self abs orption; they clas - 
sify this source as an infall candidate. ISohn et "ail (|2007f ) 
measure all three HCN (1-0) components to have two 
strong peaks, all with the blue peak higher than the red 
peak, and consider this source to be an infall candidate. 
These results are consistent with our conclusion that this 
core is likely to be approaching collapse. 
CB130: In the case of CB130, the 24 /im shadow mass 
is only marginally unstable. A plausible level of turbu- 
lent and/or magnetic support would halt collapse; for 
example, a magnetic field of ~ 70 /xG will suffice to halt 
collapse. There are no molecular line observations at the 
location of the shadow. The source CB130-3, which lies 
1 f 5 east of the shadow center, is the closest obj e ct an d 
was studied bv lPark et all (|2004l) and ILee et all (|2004fl . 
Inspection of the DSS red plate reveals what appears to 
be a cusp, or knot, of material at the CB130-3 position. 
However, our 24 /xm maps indicate that a denser and 
compact region of the CB130 complex is revealed by the 
24 /xm shadow which is completely enshrouded at near-IR 
and sho rter wavelengths (see Figure ITU]) . At the CB130-3 
position [Lee et al.1 (|2004f) measure a CS (3-2) line profile 
with a red asymmetry. Without further molecular line 
observations, we conclude that this source is probably 
stable. 

CB180: Our analysis indicates that this source is unsta- 
ble, with a 24 /xm shad ow mass that exceed s the Jeans 
mass by a factor of - 2.6. iCaselli et all (|2002f) measure an 
N2H + line width of 0.9 km s _1 but do not show the spec- 
trum of this object; this width indicates hat the contri- 
bution of turbulence is abou t 4 times that of the thermal 
support (|Hotzel et al.1 120 02:). Given the 24 /xm shadow 
mass, a magnetic field of order ~ 70 /xG would suffice 
to halt collapse. Based on these pieces of evidence, and 
without kinematical information from detailed line pro- 
files, we consider it likely that CB180 is close to collapse. 
L694-2: Our 24 /xm mass indicates that this core is sta- 
ble; however, as with L1544, L1552, and L492, this source 
has a 70 /xm shadow. Correcting our 24 /xm mass by a 
factor 5, as discussed above, yields a core mass that ex- 
ceeds the Jeans mass by a factor of ~4. Like L I 544, this 
source has been studied in some detail (e.g., ILee et all 
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the N 2 H+(2-l) and N 2 D+(2-l) line widths to be 0.27 and 
0.24 km s _1 , respectively; the N 2 H+(3-2) and N 2 D+(3- 
2) were 0.31 and 0.24 km s _s , respectively. They also 
find that this core has a large N 2 D + /N 2 H + ratio equal 
to 0.26. Similar to L1544, this core is classified by 
ICrapsi et al.l (|2005l) as highly evolve d, both chemically 
and kinematically. ISohn et alJ (2007) measure HCO (1- 
0) double-peaked l ine profiles with a blue assymctry. 
ICaselli et al.l (12001 measure a large amount of ortho- 
H 2 D + (li.o _ li,i) emission, one of the highest in their 
sample, along with L1544 (and L429). The line widths 
indicate that the turbulent support is negligible, con- 
tributing of order ~ 0.4 x the thermal support. These 
results are consistent with our conclusion that this core 
is approaching collapse. 

Here we sum marize results from two prev ious studies 
of cores CB190 (jStutz et al.ll2007t) and L429 (jStutz et all 
120091) which use similar techniques applied to 24 /xm and 
70 /xm s hadows, resp e ctivel y: 

CB190: IStutz et alJ (|2007[ ) cary out a detailed study of 
this core. The 24 /xm shadows analysis indicates that 
the mass of the core is twice the Jeans mass; however, 
the magnetic field required to support the core mate- 
rial is about 90 /xG, a plausible field strength for starless 
cores. If the magnetic field is in fact present in CB190 it 
will leak out of the core in ~3 Myr, a n ambipolar diffu- 
sion time-scale (jStahler fc P alla 2005). Turbulent pres- 
sure does not provide sufficient support against collapse. 
Millimeter line emission from various species is too faint 
to constrain the kinematics in the core. Without more 
line profile information, we conclude that CB190 is near 
collapse; 

L429: IStutz et al.1 (|2009h study the 70 /im shadow cast 
by this core and conclude that a plausible level of mag- 
netic support will not halt collapse. Additionally, the 
molecular line widths are narrow, indicating that the tur- 
bulent support is of or der the thermal support in this 
core. IStutz et al.l (I2009D conclu de that L429 is in a near 
collapse state. ILee et al.1 (|2004l ) observe a blue asym me- 
try in the CS (3 -2) line, while ICrapsi et all (|2005l ) and 
ILee et al.1 (ll9M ) measure symmetrical dou ble peaked 
N 2 H+ (1-0) and CS (2-1) lines, respectively. ISohn et all 
(2007) measure HCN (1-0) double peaked lines but note 
that this source is one of two out of 85 cores in their 
study with anomalous hyperfine line ratios. More high- 
resoluti on line observations are needed to clarify these 
results. ICaselli et all (|2008[) measure a large amount of 
ortho-H 2 D + (li,o~li,i) emission for this source, and find 
that this emission is correlated with a large degree of cen- 
tral concentration, and a large amount of CO depletion 
and deuteration, indicating that this source is evolved. 
Taken together, these pieces of evidence indicate L429 is 
close to collapse. 

In this paper we emphasize how the 24 /xm shadow ob- 
servations can reveal structure not evident in previous 
studies. We highlight the 24 /xm observations of CB130 
and L492. In the case of L492, we observe what appear 
to be 2 sub-clumps in the shadow feature (Figure [5]), ori- 
ented in a roughly east-west direction separated by about 
40", reported here for the first time. Molecular line ob- 
servations are centered on the western component (e.g., 
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iCrapsi et~a l. 2005) and indicate that this core is likely to 
be collapsing. The two components may be an example 
of core fragmentation; high spatial resolution observa- 
tions of this core should be able to confirm if both sub- 
clumps are collapsing, and to measure the relative radial 
velocity between them. In the case of CB130, we have 
discovered a 24 fan shadow which is ~ 1'.5 west of the re- 
gion prev i ously selected for millim eter line studies (e.g., 
iLee et al.l 12001 iPark et alj I2004D ; further observations 
taken at the location of the 24 /an shadow will confirm 
if this core is near collapse. 

In the top panel of Figure [55] we compare the 
24 fan shadow masses (M£f= in Table 4) to the CO- 
derived core masses (M™q c in Table 4). In all cases, 
Fcore =Mg re /Mgg e > 1.0; this discrepancy is likely due 
to the effect of freeze-out of CO in den se regions (e.g., 
iBacmann etaLll2002t [Tafalla et al.ll2002ft ; the effect will 
be stronger that illustrated for cores with 70 fim shadows 
(L1544, L1552, L492, and L694-2, see above discussion) 
because the 24 /im mass is likely to be underestimated. 
In the lower panel of FigureJ^FJwe show the 24 /jm masses 
as a function of the total molecular CO derived masses 
(integrated over a 4' radius) . Cores with 70 fim shadows 
are marked with grey arro ws indicating a pla usible mass 
correction of a factor of 5 (|Stutz et al.ll2009T ). We tenta- 
tively detect a rough trend in core mass with increasing 
envelope mass. We caution that this trend may be driven 
by distance uncertainties, small sample size, and incom- 
plete mass census. Observations from upcoming facili- 
ties, e.g., SCUBA2, which can image the sub-millimeter 
dust continuum emission, and Herschel, which can im- 
age the 70 fim shadows with a high degree of spatial 
resolution (~6") will provide further constraints on the 
evolution of starless cores. 

5. CONCLUSIONS 

We study the 8 fan, 24 /zm, and 70 /urn shadows cast 
by a sample of 14 starless cores. We derive 24 fan core 
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TABLE 1 

Spitzer IRAC and MIPS data parameters 



Source 


Other names 


R.A. 

h m s 


Dccl. 

O 1 II 




I 

O 


b 

O 


MIPS AOR a 


IRAC AOR a 


(1) 


(2) 


(3) 


(4) 




(5) 


(6) 


(7) 


(8) 


TMC1-1C 


B220 


04 41 51.4 


+25 50 


13 


174.3 


-13.3 


18155008 


5085952 


V_ LJ z o 




HA A1 97 Q 
U^t '-to z / .y 


I on AO 
-f-zy <±z 


ZD 


171^ 
± / ±.0 


in 
— 1U.D 


IZUZDoOC 


4Q1 AQAA 


CB26 


L1439 


04 59 59.8 


+52 03 


29 


156.1 


+5.9 


12020480 


4916224 


CB27 


L1512 


05 04 05.3 


+32 46 


03 


171.8 


-5.2 


12021248 


4916736 


L1544 




05 04 17.4 


+25 08 


23 


178.0 


-9.8 


18155776 


14610688 


L1552 




05 17 38.2 


+26 02 


53 


179.0 


-6.8 


18156288 


4917760 


CB42 




06 02 47.5 


+16 43 


00 


192.5 


-2.8 


12024064 




CB44 


B227.L1570 


06 07 31.0 


+19 31 


00 


190.6 


-0.4 


12022272 


4919296 


L492 


CB128 


18 15 49.1 


-03 43 


05 


25.5 


+6.1 


18159360 


14603264 


CB130 


L507 


18 16 14.3 


-02 30 


16 


26.7 


+6.6 


18159616 


5146368 


CB180 


B133,L531 


19 06 15.1 


-06 55 


25 


28.4 


-6.4 


12023808 


4925696 


L694-2 


B143,CB200 


19 41 07.4 


+10 58 


59 


48.4 


-5.8 


18160128 


14604288 



Note. — Coordinates indicate the reference pixel coordinate for each listed observation. 
AOR = Astronomical Observation Request 



TABLE 2 
HHT OTF Mapping parameters 



Core 


R.A. 


Decl. 


Obs. Date 


Size 


Line 


Tsys a RMS 6 


Vel. c 




h m s 


o / // 




' x ' 




K K-T* 


km s — 1 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) (8) 


(9) 



TMC1 


04 


41 


37.5 


+26 


02 


34 


2009 


Feb 01 


10 


X 


10 


12 co 


(2-1) 


203 





15 


2.0, 10.0 


























13 co 


(2-1) 


222 





11 


4.0, 7.0 


CB23 


04 


43 


19.1 


+29 


43 


43 


2008 


Feb 07 


13 


X 


13 


12 co 


(2-1) 


241 





22 


4.0, 11.0 


























13 co 


(2-1) 


221 





15 


5.0, 8.0 


CB26 


05 


00 


19.9 


+52 


05 


56 


2009 


Mar 02 


14 


X 


12 


12 co 


(2-1) 


263 





17 


1.7, 7.5 


























13 co 


(2-1) 


255 





25 


3.0, 7.0 


CB27B 


05 


04 


10.5 


+32 


43 


00 


2007 


Jan 28 


15 


X 


15 


12 co 


(2-1) 


218 





24 


6.0, 8.3 


























13 co 


(2-1) 


231 





25 


6.0, 8.0 


L1544 


05 


04 


16.9 


+25 


11 


50 


2009 


Feb 01 


10 


X 


10 


12 co 


(2-1) 


201 





12 


7.0, 10.0 


























13 co 


(2-1) 


220 





10 


6.0, 8.0 


L1552-SH 


05 


17 


38.6 


+26 


05 


29 


2009 


Feb 02 


10 


X 


10 


12 co 


(2-1) 


190 





15 


5.0, 10.0 


























13 co 


(2-1) 


207 





13 


7.0, 9.0 


CB42 


06 


02 


34.2 


+16 


57 


04 


2008 


Apr 16 


10 


X 


10 


12 co 


(2-1) 


264 





19 


-2.0, 9.0 


























13 co 


(2-1) 


249 





18 


1.0, 4.0 


CB44 


06 


07 


29.2 


+19 


27 


33 


2007 


Jan 28 


10 


X 


10 


12 co 


(2-1) 


193 





24 


-5.0, 3.0 


























13 co 


(2-1) 


207 





26 


-2.0, 1.0 


L492 


18 


15 


49.2 


-03 


44 


11 


2009 


Mar 01 


10 


X 


7 


12 co 


(2-1) 


299 





24 


0.0, 10.0 


























13 co 


(2-1) 


289 





33 


6.0, 9.0 


CB130-1 


18 


15 


59.8 


-02 


13 


28 


2009 


Mar 01 


12 


X 


12 


12 co 


(2-1) 


261 





19 


1.0, 11.0 


























13 co 


(2-1) 


246 





27 


5.0, 9.0 


CB180 


19 


06 


14.8 


-06 


48 


18 


2008 


Apr 18 


10 


X 


7 


12 co 


(2-1) 


209 





15 


9.0, 14.0 


























13 co 


(2-1) 


196 





15 


11.0, 13.0 


L694-2 


19 


41 


05.3 


+10 


54 


54 


2009 


Mar 01 


10 


X 


10 


12 co 


(2-1) 


297 





22 


7.0, 11.0 


























13 co 


(2-1) 


295 





18 


8.5, 10.5 



Note. — The beam size (FWHM) for the HHT at the observed frequencies is ~32" FWHM; the velocity 
resolution for these observations is 0.34 km s _1 . Objects observed in 2008 and later were observed using the 
ALMA prototype mixer, see text for details. 
a Single side-band system temperature. 
b Final OTF map RMS per pixel per velocity channel. 
c Velocity range of CO emission. 
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TABLE 3 

Shadow parameters 



Source 


R.A. 


Decl. 




Dist. a 


r b 

JVC, 24 


/(),24 


/dc,8 c 


/o,8 




h m s 


o / // 




pc 


/^Jy/arcsec 2 


/xjy/arcsec 2 


MJy/sr 


MJy/sr 


(1) 


(2) 


(3) 




(4) 


(5) 


(6) 


(7) 


(8) 


TMC1-1C 


04 41 39.1 


+26 00 


19 


140 


1390.63 


11.49 


957.92 


38.98 


CB23 


04 43 29.5 


+29 39 


11 


140 


-4.51 


3.15 


1578.37 


1.96 


CB26 


05 00 12.6 


+52 06 


02 


140 


-5.60 


2.56 


1741.00 


14.63 


CB27 


05 04 08.6 


+32 43 


23 


140 


-7.67 


5.82 


903.63 


3.88 


L1544 


05 04 16.7 


+25 10 


12 


140 


1555.05 


12.36 


4294.04 


66.20 


L1552 


05 17 38.8 


+26 05 


05 


140 


1530.99 


2.20 


1456.65 


17.89 


CB42 


06 02 27.7 


+16 56 


59 


200 


-11.0 


10.45 






CB44 


06 07 29.0 


+19 27 


42 


400 


-7.30 


3.89 


1089.62 


28.00 


L492 


18 15 48.6 


-03 45 


41 


225 


1075.20 


7.79 


2402.96 


26.16 


CB130 


18 16 11.1 


-02 16 


38 


225 


788.42 


7.30 


835.85 


28.29 


CB180 


19 06 08.4 


-06 52 


56 


400 


-5.78 


4.60 


865.37 


30.85 


L694-2 


19 41 04.1 


+10 57 


10 


230 


627.72 


6.22 


1511.84 


45.87 



Note. — The listed values of /o,8 an d /o,24 are the /rjc^subtracted quantities. Coordinates 
indicate the approximate centers of the shadow features shown in Figures [TH12I We also note 
that 1 /jjy/arcsec 2 = 23.5 MJy/sr 

a A discussion of the adopted distances and associated references is included in § 3.2. 
b /dc 24 indicated the 0.5 percentile pixel value in the 24 /an image. 
c /dc 8 indicated the 0.1 percentile pixel value in the 8 fim image. 



TABLE 4 

Core masses and sizes 



Source 


M corc 


Mggc 


a 


M co ,4 b 


p>corc c 
eff 


Dcorc c 
off 


pcorc d 






M 


M Q 




M 


arcmin 


pc 




M Q 


(1) 


(2) 


(3) 




(4) 


(5) 


(6) 


(7) 


(8) 


TMC1-1C 


1.6 




0.9 


7.1 


1.4 


0.06 


1.8 


2.7 


CB23 


2.5 




0.7 


4.9 


1.3 


0.05 


3.8 


2.0 


CB26 


5.0 




0.8 


5.0 


1.4 


0.06 


6.0 


1.6 


CB27 


2.1 




0.8 


5.6 


1.1 


0.06 


2.5 


2.4 


L15449 


1.4 




0.8 


6.1 


1.4 


0.06 


1.7 


2.8 


L15529 


8.7 




0.8 


6.8 


1.3 


0.05 


10.6 


1.1 


CB42 


0.3 




0.3 


2.9 


1.2 


0.04 


1.0 


3.0 


CB44 


35.0 




7.0 


46.2 


1.1 


0.16 


5.0 


2.8 


L4929 


3.8 




2.0 


16.1 


1.3 


0.09 


1.9 


3.3 


CB130 


3.2 




1.5 


15.6 


1.2 


0.08 


2.1 


3.1 


CB180 


10.2 




5.3 


51.5 


1.2 


0.13 


1.9 


3.9 


L694-29 


2.6 




1.7 


16.3 


1.2 


0.08 


1.6 


3.6 



a CO mass integrated over RJjg- rc , centered on the shadow coordinates. 

b CO mass integrated over an 4'-radius aperture, centered on the shadow coordinates. 

c Effective radius of the 24 fim shadow, see § 4.2.1. 

d Ratio of M=° re over Mgg°. 

c Jeans mass assuming all cores are at 10 K, using the Mj^ and R^ c to calculate the 
mean density in each core. 

f Rotational velocity required to support a core mass of M^™ for a core size of R^ rc 
s Although the quality of the 70 /im data does not allow for robust measurements, 
these cores show strong indications of having 70 (im shadows. The presence of a 70 (im 
shadow implies that the cores are optically thick at shorter wavelengths; hence M^ 10 
should be regarded as a lower limit. 
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TABLE 5 
Stability analysis summary 



G 

oOUXCC 




lyrcore l-K/rcore 
1V1 2 4 /Nlj 


/ U fiva oli&dow 


B a 


£>turb/ *^th 


V rot 

km s~ 1 


V_jI1(JII1. JZjVUI. 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


TMC1-1C 


no 


0.6 


no 








0.4 


no 


CB23 


yes 


1.3 




123 


0.3 


0.6 




CB26 


yes 


3.2 


no 
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Note. — See SS 4.3 for full discussion. 
a Magnetic field strength required to halt collapse. 

b Calculated from line widths collected from the literature, see § 4.3 for references. 
c Rotational velocity required for stability. 




Fig. 1. — 8' X 8' images at the indicated wavelengths, shown on a log scale. North is up, east is to the left. For full resolution figures 
contact stutz@mpia.de. 




Fig. 2. — Same as FigureQ] 
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Fig. 4. — Same as FigureQ] 
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Fig. 5. — Same as FigureQ] 




Fig. 6. — Same as FigureQ] 




Fig. 7. — Same as Figure[T] 
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Fig. 8. — Same as FigureQ] 




Fig. 9. — Same as FigureQ] 




Fig. 10. — Same as Figure [T] 




Fig. 12. — Same as Figure [T] 



26 





0.0 0.4 0.8 1.2 1.6 0.0 1.3 2.6 4.0 5.3 6.6 0.0 0.6 1.2 1.8 2.4 3.0 

t 24 ,2 CO T; (K km s"') 13 C0 T; (K km s"') 




0.0 0.3 0.6 0.9 1.3 0.6 1.3 2.0 2.7 3.5 4.2 0.0 0.5 0.9 1.4 1.8 2.3 

t 24 12 CO T; (K km s"') 13 C0 T; (K km s"') 



Fig. 13. — T24-maps and 12 CO and 13 CO integrated intensity maps, over the velocity ranges listed in Table 2. The black circle (1^5 
radius) indicates the area used to calculate the 24 fim shadow masses; the white circle indicates that used to calculate the molecular masses. 
The 13 CO contours levels are 0.4,0.6,0.8 times the maximum value indicated int the 13 CO color-bar. For full resolution figures contact 
stutz@mpia.de. 
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Fig. 16. — Left Column, top panel — Black histogram: 24 fim pixel value distribution for background subtracted image region of 
ro = 240" radius (except for CB42), ro is indicated in the right column as the dashed line, and is centered on the shadow. Grey histogram: 
24 /im pixel value distribution for 30" region centered on the shadow. Solid lines: minimum and maximum black histogram bins used to fit 
for the minimum a. The dashed line indicates the best-fit mean, or fo, value at the minimum best— fit value for a, indicated by the dotted 
line. Left Column, bottom panel — Best— fit values of a plotted against the maximum bin value included in the fit. The solid line indicates 
the minimum value of a. Best-fit Gaussian parameters a and the mean (fo) are indicated. Right Column: Best-fit fo and minimum a as 
a function of aperture size (ro). 
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Fig. 17. — Same as Figure IT61 
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Fig. 19. — Same as Figure IT61 



33 




20 
pixel value (MJy/sr) 




20 
pixel value (MJy/sr) 



40 



30 



20 - 



10 



x f 
+ a 



-X--X x X 



TMC1 1 C 

x x : 



+ 



+ + + + + + + + + | + + +" 



50 



100 



150 200 
radius (") 



250 



300 




150 200 
radius (") 



15 



10 



5 - 



- Xf 


i 




1 1 1 




1 




■ 


I 
1 
1 


■ ■ 1 ■ ■ 

CB26 


- + a 
















1 

1 
1 

1 


..X....X X- 










X 


X 


X 


X" 




-+ + 


+ 


+ 


X 

* + 


+ 


+ 


+ 


+ 


1 

+ 

1 
1 


+ + +: 




X 


X 












1 

1 
1 
1 




X 
















1 

1 
1 




"x . . . 


1 




1 




1 




■ 


1 





50 



100 



150 200 
radius (") 



250 



300 



Fig. 20. — Left Column, top panel — Black histogram: 8 /im pixel value distribution for background subtracted image region of 
ro = 240" radius (except for CB42), ro is indicated in the right column as the dashed line, and is centered on the shadow. Grey histogram: 
8 (im pixel value distribution for 30" region centered on the shadow. Solid lines: minimum and maximum black histogram bins used to fit 
for the minimum <r. The dashed line indicates the best-fit mean, or /o, value at the minimum best— fit value for a, indicated by the dotted 
line. Left Column, bottom panel — Best— fit values of a plotted against the maximum bin value included in the fit. The solid line indicates 
the minimum value of a. Best— fit Gaussian parameters a and the mean (/o) are indicated. Right Column — Best-fit /rj and minimum a 
as a function of aperture size ( ro). 
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Fig. 21.— Same as Figure l20l 
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Fig. 22.— Same as Figure l20l 
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Fig. 23.— Same as Figure l20l 
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Fig. 24. — Left Column — The grey scale image is 13 CO (2-1) integrated intensity map. Contours are the CO-model derived N(H) 
column density maps, which show good spatial agreement between the molecular map; contour levels are scaled to the maximum value of 
the column density in the region of the shadow. Fractional contour levels for TMC-1C are {0.8,0.9,0.95,0.995} X 1.19 X 10 22 cm -2 , for 
CB23 arc {0.6,0.8} X 1.09 X 10 22 cm" 2 , and for CB26 arc {0.5,0.6} X 1.49 X 10 22 cm" 2 . The shadow is indicated by the small white 
circle, with a radius of 1 '5; the large circle (r = 4.'0) indicates the region in which we calculate mass Mco,4 ( see Table 4 and text). Right 
Column — Mean 12 CO and 13 CO spectra for the regions indicated by the circles in the left panel. For full resolution figures contact 
stutz@mpia.de. 
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Fig. 28. — Masses presented Table 4 are plotted here. Top panel — The CO core masses calculated over the area of the core are plotted 
versus the 24 (im shadow masses. Bottom panel — The total masses (Mco,4 + MgJ 10 — Mgg ) arc plotted versus the 24 fim shadow 
masses. The arrow indicates the effect of 50% change in the assumed distance. In both panels the black curve indicates equality between 
the ordinate and the abscissa. 



